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SUMMARY 


A  Missile  performing  a  tumbling  turn  is  simulated  by  calculating  the 
angular  acceleration  produced  by  thrust  and  aerodynaaic  forces  and  inte¬ 
grating  this  acceleration  to  obtain  the  turning  rate.  The  turning  rate 
obtained  is  used  to  calculate  the  M#  matrix,  which  reflects  the  turning 
maneuver  in  the  earth  centered  inertial  (u,  v,  w)  coordinate  system. 

The  velocity  vector  turn  angle  (©)  is  defined  and  examples  are  given 
which  show  the  effects  of  initial  velocity  and  aerodynamic  force  upon 
this  angle. 

The  digital  simulation  techniques  employed  in  the  program  are  dis-  4 
cussed  as  are  the  various  program  options  and  parameters. 
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INTRODUCTION 


For  range  safety  purposes,  it  is  necessary  to  know  the  capability 
of  an  Atlas  missile  to  change  the  direction  of  its  velocity  vector.  The 
mechanism  of  this  change  could  be  autopilot  or  control  system  failure, 
the  result  being  that  the  engines  are  locked  in  some  arbitrary  fashion 
so  as  to  produce  a  turn.  The  flight  characteristics  of  an  Atlas  are 
such  that  a  condition  like  this  will  cause  the  missile  to  tumble.  The 
digital  simulation  discussed  in  this  report  is  designed  to  execute  the 
maneuvers  resulting  from  a  swivel  angle  held  constant  throughout  a 
tumbling  turn.  The  purpose  of  this  report  is  to  document  the  simulation 
and  provide  information  required  for  its  use. //The  examples  given  in 
Appendix  II  are  typical  only  and  should  not  be  considered  as  representing 
any  particular  Atlas  configuration. 
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DISCUSSION 

In  the  tumbling  turn  simulation  pr*-  .am,  it  was  necessary  to  make 
the  following  assumptions.  The  missile  is  a  rigid  body  and  can  perform 
the  maneuver.  There  is  no  fuel  or  lox  sloshing  causing  changes  in 
engine  performance  or  causing  changes  in  the  center  of  gravity.  Due  ' 
to  lack  of  aerodynamic  test  data  for  angles  of  ,a.tlack.greatcr_  than  9,0,9 1 
It  was,  necessary  to  assume  that  the  missile  has  symmetric  drag  curves; 
that  is,  the  normal  force  magnitude  would  be  the  same  whether  the  angle 
of  attack  was  10°  or  170°  and  similarly  for  the  yaw  or  side  force.  The  7 
missile  axial  drag  curve  was  modified  by  multiplying  it  by  the  cosine 
of  the  total  angle  of  attack.  At  the  start  of  the  turn,  the  missile 
is  aligned  along  the  wolocity  vector  (the  angle  of  attack  is  set  at 
zero).  A  constraint  built  into  the  program  to  satisfy  Air  Force  require- 
ments  is  that  the  tumbling  turn  be  restricted  to  a  plane.  The  program 
is  designed  to  execute  turns  in  either  the  pitch  or  yaw  plane. 

The  torques  acting  on  the  missile  to  turn  it  in  the  pitch  plane 


The  torque  due  to  the  engine  swivsl  angle 

tb  *in  6p  (xh  -  V 

,  The  torque  due  to  the  vertical  drift  of  the  center  of  gravity 
from  the  longitudinal  axis  of  the  missile 

(TB  cos  bp  +  Ts  +  Ty)  ZCG 
The  torque  due  to  aerodynamic  forces 


N  (XCG  ** 


The  corresponding  torques  for  a  yaw  turn  are: 

1.  Tb  .in  6y  <Xh  -  XCG>  =  n,  1 

2.  (T0  cos  6y  +  Tg  +  Tv)1ycg  =_  fl/J^ 

3.  y  (xCG  -  x^ .  rm  ^  A7, 
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The  torque  due  to  missile  drag  and  center  of  gravity  offpefr  is  neglected 
because  of  its  very  small  effect. 

To  find  the  angular  acf»1  oration  about  the  q  £pitch)  axis,  we  sun 
the  torques,  divide  by  the  moment  of  inertia  and  multiply  by  the  missile 
characteristic  length. 


p*  =  y—  ^  torques 


r  [».  Sln  6p(xh  -  *CG>  *  (TB  C“‘  6p  *  Ts  *  TV)  Z 

-v] 


n 

+  N  (X 


Similarly,  to  find  the  angular  acceleration  about  the  C  (yaw)  axis,  we 
may  write,  .  _/ 

y  y  /  /  ✓ 


1  = 


/  r  v  ^  y  / 

6..  -  X__)  +  (Tn  cos  6. 


/ 


[td  sin  Uy  w*h  -  ~CG'  -  '  *B - y 


y 


/ 


+  Y  <*cg  -  Xy) 


+  Ts  +  V  ycg  1 . 

J 


If  the  missile  is  constrained  to  turn  only  in  the  pitch  plane,  a 
positive  angle  of  attack  a  may  be  defined  from  the  relation, 


T_  x  1.,  *  sin  a  1 

5  v  r\ 


Similarly  for  a  turn  in  the  yaw  plane,  a  positive  angle  of  attack  3 
may  be  defined  from  the  relation, 

Ty  x  =  sin  3 

where  ly  is  a  unit  vector  in  the  direction  of  the  velocity  vector. 

Positive  angles  of  attack  will  give  positive  aerodynamic  forces. 

A  positive  pitch  swivel  angle  is  defined  to  be  one  which  will  cause 
a  turning  rate  about  the  negative  q  axis  generating  a  positive  pitch 
angle  of  attack.  A  positive  yaw  swivel  angle  6  will  give  a  turning 

y 

rate  about  the  positive  C  axis  generating  a  positive  yaw  angle  of  attack. 
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With  this  sign  convention,  we  have  the  result  that  positive  6  will 

•  •  i  P 

generate  positive  p  and  p  ;  positive  6  will  generate  positive  'i  and 

.  y 

To  obtain  the  missile  turning  rate,  we  merely  have  to  numerically 
integrate  the  angular  acceleration.  The  Adams-Bashforth  method  was  used 
to  predict  the  turning  rate  during  the  next  machine  cycle  given  three  / 

values  of  angular  acceleration. 


M  +  h  (23  i  -  16  1  ,  ♦  5  'l  J 

n  12  n  1"1  n-2 

3  1  1  1  ° 

This  turning  rarte  (with  appropriate  sign)  is  stored  for  use  in  generating 
the  M  matrix  during  the  next  machihe  cycle.  A  correction  scheme  to  obtain 


given  ^n+1  is  impractical  due  to  limitations  in  the  logical  struc¬ 
ture  of  Combo.  Accuracy  can  be  obtained  only  by  using  small  step  3izes, 
especially  in  regions  of  large  aerodynamic  forces^ 


When  the  missile  is  in  a  simulated  turn,  it  is  necessary  to  know  the 

total  angle  of  attack.  The  total  angle  of  attack  available  from  Combo 

(called  a'),  is  a  positive  angle  less  than  or  equal  to  90°.’  During  a 

tumbling  turn  however,  total  angles  of  attack  up  to  180°  are  obtained. 

During  a  360°  tumble,  a'  would  monitonically  increase  to  90®,  decrease 

to  zero,  increase  to  90°  and  decrease  again  to  zero  degrees.  What  is 

desired  is  a  function  which  would  increase  to  180°  and  th£n  decrease  to 

zero  degrees.  Such  a  function  is  derived  from  a1  and  labeled  a'  .  The 
>  P 

missile  (axial)  drag  is  multiplied  by  the  cosine  of  g1^  to  obtain  an 

approximation^ of  the  true  drag.  This  approximation  is  discussed  more 

completely  in  Appendix  III. 


The  main  purpose  behind  this  simulation  program  is  to  generate  the 
velocity  vector  turn  angle  9.  The  velocity  vector  at  the  start  of  the 
turn  is  saved  as  the  reference  vector.  During  the  turn  the  subsequent 
velocity  vectors  afe  compared  with  the  initial  velocity  vector  and  9  is 
obtained  from  the  formula, 


9 


5 
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where 

V  T *  velocity  vector  at  the  start  of  the  turn 
KjLQ 

a  any  subsequent  velocity  vector 


One  of  the  f f *  al.th.tha  ahgyy  formula -is  _that  the  aggie,  ft 

obtained  is  always  positive.  To  verify  that  such  a  condition  should  not 
always  be  true,  consider  the  following  simplified  sketch. 


Here  we  have 


•  • 


Tb  .in  6y  (Xh  -  XC(J), 


At  the  start  of  the  turn,  the  5  axis  is  aligned  along  VRIQ  and  6y  jumps 
from  zero  to  a  snail  positive  angle.  The  action  of  the  thrust  term  TR 
sin  Qy  is  to  push  the  missile  in  the  minus  q  direction,  and  add  ti  7^T, 
a  component  which  is  in  the  minus  q  direction.  The  action  of 

turning  the  missile  counter  clockwise  because  of  the  presence  of  angular 
acceleration  is  such  that  the  thrust  terra  (Tr  cos  6y  +  Tg  +  Ty)  adds  to 
V  a  component  A  V  T  in  the  positive  q  direction.  As  it  turns  out,  the 
velocity  vector  turns  first  clockwise  and  then  counter  clockwise  in  the 
above  sketch.  This  is  due  tc  the  lag  of  the  missile  in  getting  started 
in  its  turn.  Since  VRJ  at  some  time  t  after  the  start  of  the  turn  crosses 
Vjjiq  t  *11  values  of  0  up  to  this  time  are  defined  to  be  negative. 
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PROGRAM  STRUCTURE 

The  tumbling  turn  program  is  divided  into  four  parts.  The  first  part 
is  a  one  shot  pass  to  set  up  initial  conditions.-  The  second  part  is  sub¬ 
routine  70  which  sets  up  the  engine  thrust  along  the  missile  axes.  Part 
three  sets  as  a  guidance  program  to  turn  the  missile.  The  final  paTt  of 
the  program  is  a  termination  code  to  end  the  run. 

In  the  initial  conditions  set  up,  various  flags  and  storage  locations 
are  zeroed.  Also,  the  matrix  is  changed  to  gero  out  the  angles  of  attadc 
q  and  0.  This  is  accomplished  by  replacing  the  matrix  with  where; 


COSCX 

0 

sina 

cosB 

sin0 

0 

0 

1 

0 

•  M0  = 

-sin0 

COS0 

0 

-sina 

0 

cosa 

0 

0 

1 

Combo  subroutine  70  sets  up  engine  thrust  along  the  missile  axes.  If 
a  pitch  turn  is  desired,  zero  thrust  is  along  the  q  axis  and  the  thrust 
along  the  £  axis  is  -Tfl  sin  6^.  If  the  turn  is  in  yaw,  -Tg  sin  6^  is  the 
thrust  along  the  tj  axis  and  zero  thrust  is  in  the  C  direction.  In  either 

case,  the  thrust  along  the  ^  »xis  is  Tg  cos  6  +  Tg  +  Ty. 

*  *  \ 

In  the  guidance  portion  of  the  tumbling  turn  program,  the  remaining 

variables,  angular  acceleration,  turning  rate,  total  angle  of  attack,  and 
the  velocity  vector  turn  angle  are  calculated.  It  is  at  this  section  of 
Combo  that  the  position  anti  velocity  of  the  simulated  missile  are  correct 
and  all  necessary  inputs  to  the  tumbling  turn  program  are  available. 

After  deciding  whether  the  turn  is  in  pitch  or  yaw,  the  torques  acting  on 
the  missile  are  calculated  and  summed.  Moment  of  inertia  is  found  from  a 
table  of  moment  of  inertia  vs.  weignt.  To  obtain  the  turning  rate,  the 
following  procedure  is  used.  If  the  cycle  number  is  one,  angular  accel¬ 
eration  is  assumed  constant  and  the  turning  rate  is  found  by  multiplying 
the  angulur  acceleration  by  the  step  size.  If  the  cycle  counter  is  two, 
the  rate  of  change  of  angular  acceleration  is  assumed  constant  and  a  simple 
two  point  formula  was  derived  to  obtain  the  turning  rate.  When  the  cycle 
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counter  is  three  or  higher,  the  Adama-Bashforth  formula  take#  over.  The 
total  '"'gle  of  attack  and  the  velocity  vector  turn  angle  are  obtained  as 

/ 

outlined  earlier. 

Termination  subroutine  19  is  designed  to  end  all  calculations  when  the 
maximum  value  of  the  velocity  vector  turn  angle  is  obtained. 
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PROGRAM-INPUTS 


A  referent*  trajectory  must  be  available  to  specify  the  position  and  ' 
velocity  vectors,  thrust,  weight,  the  M(  matrix  and  time  as  inputs  to 
Combo  at  the^time  a  tumbling  turn  is  to  be  a tarted.y^A~*table  of  moment  of 
inertia  vs,  weight  must  bo  entered  also.  At  the  present  time,  room  for 
only  on*  table  of  40  entries  is,  available ;  therefore,  il  the  pitch  moment 
of  inertia  differs  a  great  deal  from  the  yaw  moment  of  inertia,  only  one 
type  of  turn  will  be  peseible  without  changing  the  table.  Aerodynamic 
functions  should  be  entered  as  tables  rather  than  polynomials. 

Four  inputs  specify  the  tumbling  turn  program  parameters.  The  engine 
swivel  angle  6  is  entened  in  ENGSA.  PYAtfF  specifies  a  pitch  or  yaw  turn. 

A  zero  in  this  location  identifies  a  pitch  turn,  any  positive  quantity 
specifies  a  yaw  turn.  INMPY  and  INADD  specify  moment  of  inertia  multiply¬ 
ing  and  additive  factors  respectively^  which  are  useful  in  making  disper¬ 
sion  runs. 

If  termination  option  19  is  jused,  a  non-zero  positive  quantity  must 
be  specified  in  STGK.  This  number  represents  the  elapsed  time  that  takes 
place  before  testing  for  a  peak  9  value.  This  is  necessary  in  order  to 
bypass  the  region  of  negative  0's  (defined  previously)  where  the  logic  of 
this  routine  breaks  down.  Also,  by  judiciously  picking  this  number,  it 
would  be  possible  to  terminate  on  the  second  or  later  peaks.  A  typical 
value  is  STGK  =  2. 

A  card  punch  output  is  included  in  the  tumbling  turn  program.  This 
is  convenient  for  both  tabular  listings  and  automatic  graph  plotting.  Two 
cards  are  punched  each  machine  cycle.  On  the  first  card,  time  and  the 
position  and  velocity  vectors  in  downrange  coordinates  ore  punched  in 
fixed  point  arithmetic.  On  the  second  card,  tba  quantities  punched  are 
time,  the  velocity  vector  turn  npgle,  .uc  v*  cvvacic,  9pucd, 

elapsed  time,  and  the  engine  swivel  angle  (pitch  and  yaw). 
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CONCLUSIONS 

The  accuracy  obtained  from  this  program  is  quite  high  when  the  turns 
take  place  out  of  the  atmosphere.  Turns  starting  in  the  atmosphere  should 
be  within  accuracy.  A  disadvantage  of  this  program  is  the  high  cost  of 
computer  time  because  it  runs  relatively  slowly  (small  integration  step 
size) . 
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APPENDIX  I 


DEFINITION  OF  SYMBOLS 


Symbol 


Definition 


Units 


Integration  stepsize 

Moment  of  Inertia  about  pitch  axis 


seconds 


slugs/ft, 


Moment  of  Inertia  about  yaw  axis 


slugs/ft.' 


Missile  characteristic  length 


125  ft. 


Matrix  which  transforms  a  vector  from 
missile  to  inertial  coordinates 


Aerodynamic  normal  force 


Booster  engine  thrust 

Sustainer  engine  thrust 

Vernier  engine  thrust 

Missile  initial  velocity  vector  with 
respect  to  wind  in  inertial  coordinates 

Missile  velocity  vector  with  respect  to 
wind  in  inertial  coordinates 


Axial  location  of  center  of  gravity  fro* 
missile  reference  line 


unit  of  L 


Engine  hinge  point  location  from  missile 
reference  line 


unit  of  L 


Normal  force  center  of  pressure  location 
from  missile  reference  line 


unit  of  L 


Side  force  center  of  pressure  location 
from  missile  reference  line 


unit  of  L 
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Symbol 


a 


o' 

fi 


6 


P 


6 


y 


o 


p 

•  » 
p 


APPENDIX  I  (Cont'd) 

Definition 

Aerodynamic  side  force 

Side  drift  location  of  center  of  gravity 
from  longitudinal  axis.  Positive  in  direc¬ 
tion  of  +i}  axis 

Vertical  drift  location  of  center  of  gravity 
from  longitudinal  axis  of  missile.  Positive 
in  direction  of  +(  axis 

Pitch  angle  of  attack  (positive  nose  up) 

Total  angle  of  attack 

Total  angle  of  attack  derived  from  o' 

(see  text) 

Yaw  angle  of  attack  (positive  nose  left) 

Turning  rate  about  yaw  axis 

Angular  acceleration  about  yaw  axis 

Engine  swivsl  angle  in  pitch  plane. 

Positive  for  nose  up 

Engine  swivw  1  angle  in  yaw  plane. 

Positive  for  nose  left 

Velocity  vector  turn  angle. 

(See  text  for  definition  of  sign  convention) 

Turning  rate  about  pitch  axis 

Angular  acceleration  about  pitch  axis 


/ 


Units 

lbs. 

unit  of  L 

unit  of  L 
degree 
degrees 
degree 

degree 

rad/sec 

rad/sec^ 

degree 

degree 

degree 

rad/sec 

A/  2 

rad/sec 
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APPENDIX  II 
EXAMPLES 


In  order  to  aid  the  prospective  user  of  this  program  to  visualise 
what  to  expect  for  results,  graphs  have  been  prepared  for  an  arbitrary 
Atlas  configuration  undergoing  a  tumbling  turn  at  three  different  times 
during  powered  flight.  The  times  selected  are  early  in  booster  phase 
(16  sec.),  midway  in  booster  with  high  dynamic  pressure  (64  sec.),  and 
late  in  sustainer  phase  (250  sec.). 


Figure  3  (a)  shows  a  graph  of  the  velocity  vector  turn  angle  0  vs. 
elapsed  time  16  seconds  after  launch  with  various  engine  swivel  angles. 
The  predominant  characteristics  are  the  large  values  of  0  and  the  non¬ 
linear  envelope  of  the  peak  values  of  0.  The  large  values  of  0  result 
because  the  initial  velocity  vector  is  relatively  small  in  magnitude  and 
the  velocity  added  during  the  turn  is  a  high  percentage  of  this  initial 
velocity.  Suppose  that  at  the  start  of  the  turn,  the  initial  velocity  is 
V  T_  and  that  the  velocity  added  during  the  turn  is  AV.  Then 


VRI0  x  ^rI6'/+ 


AV)| 


sin  9^  * 


RIO 


VRI0  +  *V 


V 

_RI0 

X 

AV 

* 

v 

RIO 

VRI0  *  *V 

Now  imagine  that  the  initial  velocity  is  increased  to  aVRIQ  ,  «  >  1»  but 
the  same  velocity,  AV,  is  added  during  the  turn.  Then, 


sin  0, 


*VRI0  X  (*VRI0  * 
Kio!  I-RIO  ♦ 


>1 


I* 


RIO 


x  AV 


RIO 


*VRI0  *  *V 
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While  it  is  net  true  that  the  velocity  gained  during  a  turn  will  be  the 
same  at  different  start  up  times,  we  see  from  this  simplified  analysis 
that  the  general  trend  is  toward  smaller  values  of  6  with  greater  values 


of  V 


RIO* 


The  smaller  the  value  of  the  swivel  angle  6,  the  greater  will 


be  the  added  component  of  velocity  AV.  It  would  appear  then  that  with  a 


given  initial  value  of  that  the  smaller  value  of  6  would  give  the 

greater  value  of  0.  A  glance  at  the  grstph  shews  this  to  be  only  partly 


true.  Designating  0^  the  maximum  value  of  0  obtained  with  a  swivel  angle 
of  5®,  0g  the  value  obtained  with  6=3°  etc.,  we  have 


85<83<82'  80<81<82 

Essentially,  it  is  possible  to  optimize  the  velocity  vector  turn  angle 
through  proper  choice  of  the  swivel  angle,  a  most  curious  result. 

This  dilema  can  be  explained  in  the  following  manner.  The  most 
significant  term  in  the  equation  to  calculate  angular  acceleration  is  the 
aerodynamic  moment,  but  only  after  enough  angle  of  attack  has  been  gone* * 
rated  by  the  moment  terms  containing  thrust  and  swivel  angle.  Once  a  small 
angle  of  attack  is  obtained,  the  aerodynamic  moaent  becomes  larger  and 
larger,  generating  larger  and  larger  angles  of  attack,  in  other  words,  the 
angle  of  attack  and  the  aerodynamic  moment  are  regenerative  (up  until  the 
angle  of  attack  is  90°).  The  smaller  the  swivel  angle,  the  longer  the 
period  of  elapsed  time  will  be  until  the  aerodynamic  moment  can  take  over 
and  control  the  turn.  During  the  elapsed  time  before  the  aerodynamic 
moment  becomes  appreciable  in  magnitude,  the  velocity  vector  is  increased 
in  magnitude  with  little  change  in  direction.  An  equivalent  turn  could  be 
obtained  by  starting  the  turn  later  in  flight  but  using  a  larger  swivel 
angle.  But  it  was  shown  that  by  increasing  vRjq  »  the  maximum  value  of 
<  9  will  be  made  smaller.  The  reasoning  envolved  in  this  discussion  can  be 
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Bade  apparent  by  referring  to  Figure  3  (b).  This  graph  shows  the  angular 
acceleration  of  the  missile  during  the  turn.  The  acceleration  curve  for 
the  zero  degree  swivel  angle  is  much  greater  than  the  others,  showing  that 
the  velocity  of  the  missile  is  much  greater  at  equal  angles  of  attack. 
Figure  3  (c)  is  a  graph  of  the  Bissile  turning  rate  with  various  swivel 
angles.  Notice  that  the  peak  values  of  turning  rate  are  inversely  propor¬ 
tional  to  the  peak  values  of  0  at  constant  swivel  angle. 

Figure  4  (a),  Figure  4  (b)  and  Figure  4  (c)  are  graphs  of  the  velocity 

vector  turn  angle,  angular  acceleration,  and  turning  rate  respectively  for 
a  turn  in  a  gegion  of  very  high  dynamic  pressure.  Notice  that  the  peak 
values  of  all  these  curves  are  almost  independent  of  engine  swivel  angle. 
Also,  the  elapsed  time  is  very  short  as  a  consequence  of  the  very  high 
aerodynamic  forces. 

Figure  5  is  a  graph  of  the  velocity  vector  turn  angle  during  sustainex* 

engine  phase.  The  only  moments  acting  to  turn  the  missile  are  due  to 

thrust,  the  aerodynamic  forces  being  zero.  The  angular  acceleration  was 
not  plotted  because  it  was  practically  constant  throughout  the  turn  for 
each  of  the  swivel  angles.  Correspondingly,  the  turning  rates  were  linear 
functions  with  slope  dependent  on  the  swivel  angle.  An  important  charac¬ 
teristic  of  the  class  of  turns  taking  olace  out  of  the  atmosphere  is  that 
the  envelope  of  the  curves  for  various  swivel  angles  at  any  particular 
start  up  time  is  a  straight  line. 
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APPENDIX  III 

TESTS  OF  ASSUMPTIONS 


As  stated  earlier  in  the  report,  certain  assumptions  were  made 
regarding  the  aerodynamic  characteristica  of  the  Atlaa  missile  because  of 
the  lack  of  data.  0*e  of  these  was  that  the  missile  axial  drag  curve  was 
modified  by  multiplying  it  by  the  cosine  of  the  total  angle  of  attack  to 
account  for  the  change  in  axial  drag  with  the  large  angles  of  attack 
encountered  during  a  turn.  Another  assumption  was  that  the  normal  and  side 
forces  on  the  missile  were  symmetric  about  the  angle  of  attack  equal  to 
90*.  To  test  these  assumptions,  others  were  made  in  their  place,  the  i^ea 
being  that  if  only  small  effects'  were  observed,  the  original  results  could 
be  considered  quite  satisfactory. 

To  test  the  original  assumption  about  drag  varying  as  the  cosine  of 
the  angle  of  uttuck,  a  trial  run  was  made  assuming  that  the  drag  should  be 
multiplied  by  the  square  of  the  cosine  of  the  angle  of  attack.  This  was 
tested  only  for  the  turn  taking  place  with  very  high  dynamic  pressure 
(64  seconds)  as  it  was  felt  that  this  would  present  the  greatest  variation 

I 

from  the  original  results.  The'  peak  value  of  9  obtained  under  these  condi¬ 
tions  occurred  at  the  same  elapsed  time  and  had  an  increase  in  amplitude  of 
.46*. 

To  teat  the  original  assumptions  about  a  symmetric  normal  force,  a 
non-symmetric  one  was  tried  and  two  tests  made,  one  at  16  seconds  and  ono 
at  64  seconds.  The  original  normal  force  (N)  was  multiplied  by  the  quantity 


It  is  easily  seen  that  in  the  region  of  uonotonic  increasing  a’,  0  a' 
90*,  the  multiplier  m  (a1,  a^')  =  *  because  a'  =  a^'.  Once  past  tjiis 
point  however,  becomes  the  supplement  of  a'  and  m  (a',  0^')  increases 

limearly  to  its  maximum  value  m  (0,  180*)  =  2.  With  the  normal  force 


multiplied  by 


),  the  following  results  were  obtained.  The  peak 
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value  of  9  obtained  from  the  16  sec.  case  was  1.5794  lew  and  occurred  .92% 
sooner  in  elapsed  time  when  compared  with  the  reference  (unmodified  N) 
case.  The  peak  value  of  9  obtained  from  the  turn  starting  at  64  sec.  was 
.8294  low  and  occurred  1.2894  sooner  as  compared  with  reference.  The  direc¬ 
tion  of  the  shift  of  9  in  both  cases  is  such  that  the  envelope  of  the  turns 
for  varisus  swivel  angles  would  be  little  changed. 
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